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ABSTRACT 

We show that large, high-redshift (z > 10) galaxies with virial temperature 
in excess of 10 4 K may be mostly comprised of cold atomic clouds which were 
formerly minihalos. These clouds move at a speed of ~ 15 — 30km s" 1 and collide 
with one another on a time scale of 10 7 yr. The supersonic collisions result in 
spatially distributed star formation with high efficiency We then show that most, 
subsequent star formation in cold atomic clouds may be triggered by Shockwaves 
launched from the first stars formed in collisions. Those shockwave-compressed 
clouds are even more widespread spatially because of large imparted velocities 
and some can possibly escape into the intergalactic medium. More importantly, 
with respect to cosmological reionization, widespread star formation would allow 
a much higher ionizing photon escape fraction. These favorable conditions may 
form the physical basis to enable the standard cosmological model to produce 
a reasonably high Thomson optical depth r e = 0.10 — 0.14. A chain reaction 
of star formation in minihalos in the intergalactic space may be triggered by 
explosions in the intergalactic medium, if minihalos are strongly clustered. In 
this still higher r e would be achievable. 

Subject headings: reionization - Shockwaves - collisions - cosmic microwave 
background 
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1. Introduction 

The Wilkinson Microwave Anisotropy Probe (WMAP) polarization observations 
measured Thomson optical depth r e = 0.17 ± 0.04 (68%) (Kogut et al. 2003), requiring 
a full reionization to have occured first at z = 20ig°, such as in the double reionization 
picture (Cen 2003a). This high Thomson optical depth is further confirmed by a variety of 
combined data sets (Tegmark et al. 2003). To ionize the universe at such a high redshift 
demands both a high formation efficient (e) for Population III stars and a large escape 
fraction (f esc ) for the produced ionizing photons. Quantitatively, one would require the 
product of the two to be greater than 0.01 in order to produce r e > 0.11 (e.g., Cen 2003b; 
Haiman & Holder 2003; Wyithe & Loeb 2003; Ciardi, Ferrara & White 2003; Ricotti & 
Ostriker 2003) to be consistent with the WMAP observations at ~ 1.5a level. Conventional 
wisdom may view these assumptions somewhat too optimistic at first sight. 

In this Letter we suggest that star formation at high redshift [z > 10) may occur 
in massive starbursts, triggered by supersonic collisions between constituent cold atomic 
clouds and blastwaves produced by explosions of first stars. The consequences are (1) high 
star formation efficiency and (2) high ionizing photon escape fraction. 

2. Collisions between Cold Atomic Gas Clouds in High Redshift Galaxies 

The slope of density power spectrum is ~ —3 in the cold dark matter cosmological 
models. As a result, galaxies over a wide range of mass, 10 4 — 10 8 M of interest here at 
high redshift, are expected to form nearly simultaneously, when they form. This would 
imply that large galaxies of mass 10 7 — 10 s M Q (virial temperature T v > 10 4 K at z > 10), 
have most of their mass in small clouds, which are formerly minihalos. How much mass 
of a large galaxy is in minihalos? This is computable approximately using the Extended 



-4- 



Press-Schechter (EPS) theory (Bond et al. 1991; Bower 1991; Lacey & Cole 1993, LC 
hereafter). We use the instantaneous merger rate for a halo of mass M p to merger with a 
halo of mass AM = M — M p to form a halo of mass M (LC): 

r-, r\ Ti j / 2 x1/9 l,dln5 ci 1 , diner . 5 c (t) 1 r 5 c (t) 2 , 1 1 

R N (M V -> M)dM = — ; 0Nq/0 exp — — ~)]dM, 

NK v ! V V dint l M l d\nM l a (1 - a 2 /a*) 3 / 2 PL 2 V 2 a 2 ' 1 ' 

(1) 

where a is the density variance of a sphere in the linear density field extrapolated to z = 
that contains mass M; a p is <j for mass M p ; S c (t) is the critical density for collapse (equation 
2.1 of LC); t is age of the universe at the epoch examined. The total mass accretion/merger 
rate of halos of mass in the range M 1 to M 2 to form a halo of mass M at t is then 

rM-Mi 

R M (M, Ml, M2) = / R N (M P -> M)(M- M p )dM p . (2) 

JM~M 2 

Figure 1 shows the mass fraction of accreted objects in the mass range 2 x 10 4 — 10 7 M by 
a large halo with M — 1 x 10 s M as a function of redshift. We see that 50 - 80% of all 
mass in large galaxies with virial temperature greater than 10 4 K at z > 10 are in minihalos. 
The standard WMAP-normalized cold dark matter model (Spergel et al. 2003) is used: 
Q M = 0.27, A = 0.73, Q b = 0.04, H = 70km/s/Mpc, n = 0.99 and a 8 = 0.90. 

The constituent minihalos do not remain intact in the process of infall and virialization. 
First consider the effect due to gravitational tidal field of the large, host halo on minihalos. 
For simplicity we will assume that all dark matter halos are singular isothermal spheres 
truncated at the respective virial radii, p(r) = p v (r/r v )~ 2 for r < r v , where p v = lOOp, p is 
the mean density of the universe at that redshift and r v is the virial radius. For illustration 
we will consider a large halo of mass 1 x 10 8 M at z = 20 having a virial temperature of 
3.7 x 10 4 K. We assume that 2/3 of the total mass of the large halo were originally from 
minihalos (see Figure 1), and all minihalos have the same mass of 1 x 1O 6 M with a 
virial temperature of 1700K, totaling 67 such minihalos in the large halo. We denote the 
virial radius of the large halo as tl (= 0.71kpc proper) and that of the minihalos as tm 
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Fig. 1. — shows the fraction of total mass acquired by a halo of mass M — 1 x 10 8 M 
at redshift z that is in minihalos in the mass range 2 x 10 4 - 1O 7 M . Note M = 
1 x 1O 8 M corresponds to a halo with a virial temperature of 3.7 x 10 4 K and halos with 
M = 2 x 10 4 - 1O 7 M have a virial temperature of 120 - 7900K, i.e., minihalos, at z = 20. 
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(= 0.15kpc proper). Then, a minihalo at radius r from the center of the large halo would 
have its matter exterior to r T = r M (r/r L ) be tidally stripped off. The mean collision rate 
between minihalos (ignoring gravitational focus due to high speed) is 

Jq L vn M (r)(r(r)n M (r)4Trr 2 dr 

< nav >= — — — , (3) 

Jo n M (r)^nr 2 dr 

where a(r) = ~kt\ is the cross section of a minihalo at radius r; um{t) is the number density 
of minihalos at radius r; v is the velocity of each minihalo. Assuming also an isothermal 
distribution of minihalos, Um{t) = n (r /r^) -2 , we have 

< nav >= 7cr 2 M vn . (4) 

This result is obtained only considering the gravitational tidal effect. Let us now consider 
the effect due to external inter-minihalo medium (IMM) on the internal gas of minihalos. 
Consider a minihalo located at r L . The IMM with T = 3.7 x 10 4 K and a density of p v /3 has 
a pressure higher than the minihalo gas with T = 1700K and a density of p v . Hence, the gas 
density in minihalo would be raised by adiabatic compression by the external gas by a factor 
of 3.3 and the temperature increases by a factor of 2.2, resulting in a pressure equilibrium 
(ignoring self-gravity of a minihalo for this simple calculation). The compression would 
reduce the cross section by a factor of 2.2 (A better treatment will be using virial theorem 
with external pressure but is not given here). The percentage reduction in cross section at 
other radii is the same. Thus, the average collision time between minihalos is 

tc = a' 2 = 43 x 10 V, (5) 

where the minihalo clouds in the large halo are assumed to move at a velocity equal to 
v = \/3a v = 30km s _1 (a v = 17.4km s _1 is 1-d the velocity dispersion of the large halo for 
our illustrated case). It is seen that the collision time among minihalos is about 25% of the 
age of the Universe at z = 20. The time it takes for a minihalo to cross the large halo is 
4.6 x 10 7 yr. Thus, on average, a minihalo experiences a collision for every crossing time. 
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Gas within minihalos may be stripped off by ram pressure when traveling through 
the IMM of the large host halo. Applying the same physical argument of Gunn & Gott 
(1972), if the gravitational restoring force of a minihalo is exceeded by ram pressure force, 
gas and dark matter of the minihalo will be displaced, i.e., gas being stripped off. This may 
be expressed by the condition Pimm v q = Pmv m , where Pimm is the density of the IMM; 
Vq is the traveling velocity of the minihalo; pu is the mean density of a gas parcel of the 
minihalo in question; % is the 3-d velocity dispersion of the minihalo. We see that this 
condition is similar to the pressure equilibrium condition above (without considering the 
self-gravity of minihalos). Thus, ram pressure stripping is relevant and likely to act at the 
same time when the cloud is compressed by external thermal pressure. It is possible that 
gas in outmost layers of a minihalo may be subject to ram pressure stripping while most 
interior gas will not be dislodged. We will ignore the change to tc (equation 5) possibly 
caused by ram pressure stripping, because other effects such as deformation of minihalo gas 
clouds may easily dominate but are hard to characterize simply. 

What would happen in a head-on collision between two clouds of cold gas? Smith 
(1980) and Struck-Marcell (1982a,b) have studied this issue of cloud collision inside galaxies 
for different ranges of Mach numbers using both analytic method and simulations, in 
somewhat different environments. Here we focus on the thermodynamic properties of 
post-shock gas, computed assuming that post-shock region is isobaric, while leaving detailed 
treatment to future 3-d hydrodynamic simulations. All cooling/heating processes for a 
primordial gas are incorporated, including atomic cooling, H 2 formation and cooling and 
Compton cooling with the rates from Abel et al. (1997). Radiation background is assumed 
to exert no significant effect, as shown by Mac Low & Shull (1986) in the context of 
intergalactic shocks, which are physically similar to collision-induced shocks considered here. 
Figure 2 shows the evolution of the postshock gas following a head-on collision between two 
cold atomic gas clouds. Note that the time for the shock to propagate through an entire 
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Fig. 2. — shows the evolution of the postshock gas during a head-on collision between two cold 
atomic gas clouds. The panels from top to bottom show the evolution of has temperature, 
density, molecular hydrogen fraction and Jeans mass, respectively. Two cloud each has an 
initial velocity of Vq = 30km s -1 in the lab frame with opposite direction with pre-collision 
parameters: T = 1700K and n = 0.42cm -3 . The vertical dashed tick sitting at the bottom 
axis of the bottom panel indicates the elapsed time of 0.5t s = 3.64 x 10 6 yr, since the onset 
of the collision. The black horizontal dotted line in the bottom panel indicates the total 
baryonic mass of the collided clouds. 
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cloud is t s = 2r M /(4f /3) = 7.3 x 10 6 yr. Since the shock heated cloud has already cooled 
to a much lower temperature and the Jeans mass has become much smaller than the total 
baryonic mass of the merged cloud at a time scale shorter than for the shock to transverse 
the entire cloud, the subsequent rarefaction wave following the thermalization (shock) 
process is insignificant, so are the compression and expansion phases following that. This 
indicates that the isobaric assumption is approximately valid and our simple calculation 
should have captured the essence of the process. We see that a molecular hydrogen fraction 
of 8 x 1CT 4 is reached and gas is cooled to ~ 100K within a fraction of the Hubble time 
after the collision (~ 10 6 — 10 7 yr). The final Jeans mass is about 100 M , indicating that 
stars formed may be massive. One interesting difference in terms of star formation between 
this picture and that of Abel et al. (2001) is that multiple stars, possibly a star cluster, may 
form, which, unlike in a single star case of Abel et al. (2001), may provide an additional 
mechanism to remove the angular momentum of the collapsing gas through orbital coupling. 

In summary, supersonic collisions between cold atomic gas clouds within large halos 
can trigger star formation and star formation efficiency may be quite high, because a large 
portion of the colliding clouds may be able to cool and form stars. The spatially distributed 
star formation will be much more favorable to the escape of ionizing photons. 

3. Star Formation Triggered by Shockwaves 

The possible spread in time among collisions may be 10 7 — 10 8 yr, longer than the 
lifetime of massive first stars of ~ 3 x 10 6 yr. Thus, it is likely that collisions between 
minihalo clouds can not be considered coeval. Consider two possible scenarios. 

In the first scenario, the first stars were so massive that they end their lives imploding 
and display no explosive events, for example, if the stars are more massive than 300 M Q and 
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non-rotating (Rakavy, Shaviv, & Zinamon 1967; Bond, Arnett, & Carr 1984; Glatzel, Fricke, 
& El Eid 1985; Woosley 1986), the subsequent star formation in collisions that follow would 
proceed normally. Star formation efficiency should be high and so is the ionizing photon 
escape fraction, since star formation will be spatially distributed and extended (see below). 

In the second scenario, first stars end in explosions and launch a blastwave. Before 
computing the dynamics of blastwaves from the explosion of first stars, let us first consider 
the photoionization effect by the first stars on non-star-forming minihalos, before blastwave 
encompassing them. The ionizing photons from the massive stars formed in collisions will 
propagate outward from the minihalo. Since the ionization front (assuming an ionizing 
photon emissivity of 1.6 x 10 48 photons s -1 M _1 ; Bromm, Kudritzki, & Loeb 2001) travels 
at a velocity of v IF = F/n = 3478(rf/lkpc)^ 2 (e/0.1)(n/0.1cm _3 ) _1 km/s (assuming no 
attenuation up to distance d from the emitting source), much faster than the sound speed 
in photoionized gas, an R-type ionization front thus develops. Subsequently, a sound wave 
(not a shockwave) will develop driving the ionized gas outward at the sound speed, which 
would travel 0.046kpc (~ rjw/3) in the lifetime of the massive stars. As a result, the gas 
density immediately surrounding the stars will be much reduced. Given that the expansion 
distance is comparable to the original minihalo radius one may expect that the final gas 
distribution surrounding the stars at the end of their lives may have a density comparable 
to mean density of the original minihalo. Without detailed simulations we resort to some 
order of magnitude estimates. Let us assume that the clumping factor of the gas in 
the star forming minihalo cloud is C (in terms of pre-shock gas density), then it can be 
shown that the total amount of gas required to balance the ionizing photon emission from 
recombination is M rec = 1.0 x 10 5 (e/0.1)(C/16) _1 M for an ionizing photon luminosity of 
4.2 x 10 52 (e/0.1)sec _1 , which should be compared to 2.6 x 1O 5 M of the total gas mass of 
the colliding clouds. The fact that the former is of the same order of the latter and the 
fact that the post-collision gas medium of the minihalos is likely to be quite irregular tell 
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us that (1) a significant fraction of the ionizing photons is likely to escape the star- forming 
minihalos, and (2) the Stromgren sphere size is not much larger the minihalo itself, which 
would imply that the ionizing effect from the star-forming minihalos on other minihalos 
in the same galaxy may be relatively mild. In a fraction of the solid angle, fa, from the 
star-forming minihalo optical depth may be small enough to allow for the ionizing photons 
to escape. For those directions, ignoring the recombination effect due to the IMM, it can be 
shown that the minimal distance from the star-forming minihalo for the recipient minihalo 
of ionizing photons to experience a D-type ionization front is d m i n = 2.2(-^-) _1/,2 kpc, which 
should be compared to the virial radius of the large-host halo of r L = 0.71kpc. Thus, within 
the fraction of the solid angle where ionizing photons are able to escape into the IMM, it 
appears that ionization front is R-type initially at the virial radius of a minihalo but turns 
into a D-type front when penetrating deeper into the minihalo (at 0.57tm for a minihalo at 
a distance of r^). A D-type front may transmit a weak shock causing slight compression of 
the gas, while the regions cruised through by the R-type front may expand outward. In 
brief, we see that a small fraction of non-star-forming minihalos within the solid angle of 
escaped ionizing photons from star-forming minihalos may experience some dynamic effect 
due to photoionization, likely resulting in a very compact central region and evaporated 
outskirts (Barkana & Loeb 1999; Shapiro et al. 2003). Most of the minihalos are, however, 
not irradiated by ionizing photons from the first stars and would remain neutral and intact. 

Now consider the effects of the blastwaves produced by the explosions of the first stars 
in the same large host galaxy. We will assume that the first stars each have mass of 150 M 
with supernova explosion energy of 3 x 10 52 erg. A star formation efficiency of e would give 
a total explosive energy of 5.2 x 10 54 (e/0.1)erg. It is not required that the mass of the first 
stars be in the range 130 — 300 M Q to induce pair-instability supernovae; our conclusions 
remain little changed if the mass of the first stars is the range of 20 — 130 M . We will 
now estimate how far and at what speed the blastwave from such an explosion would 
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travel. For simplicity and easy tractability we assume that the blastwave has two phases: 
the Sedov- Taylor phase followed by a momentum conserving phase (with an instantaneous 
cooling phase in-between). The assumed second phase is conservative in the sense that it 
ignores any pressure of possibly still hot interior gas on the shock shell. The formula for the 
first phase is adopted from Ostriker & McKee (1988) assuming a uniform density locally. 
Cooling rate for a primordial gas is taken from Sutherland & Dopita (1993) with Compton 
cooling added. Figure 3 shows the velocity of the blastwave after it has traveled a distance 
such that it would, on average, hit 1 (solid curve) and 8 (dashed curve) other minihalos, as 
a function of the ambient gas density (depending upon the radius). We see that for density 
in the relevant range the expected shock velocity sweeping minihalos may be in the range 
v s = 20 — lOOOkm/s. We note that, with significantly reduced density surrounding the 
exploding stars due to photoionization, the blastwave only starts cooling way after breaking 
out the minihalo gas clouds in which the stars were formed. Thus, neglecting the gas in the 
host minihalo does not introduce significant errors. 

What would happen to the minihalos when shocks of this magnitude travel through 
them? A shock with a velocity of 93km/s in the ambient medium would produce a shock 
of velocity 30km/s (= 93 (pimm / Pm) 1 ^ 2 = 93/v / 9~l3) upon entering a minihalo, producing 
results comparable to those shown in Figure 2. One important difference is that, in this 
case of encompassing blastwaves, the external thermal pressure is high and could aid in 
compressing the shocked and cooled minihalo gas clouds, which will likely enhance star 
formation in the minihalos. The reader may consult Kang et al. (1990) for a 1-d treatment 
in the context of globular cluster formation and Klein, McKee, & Colella (1994) for a recent 
calculation. In short, blastwaves from first massive stars could cause star formation in 
other minihalos in the same galaxy, yielding a widespread star formation. Most likely, one 
generation of stars may be able to trigger star formation in majority of the minihalos. 
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Fig. 3. — shows the velocity of the blastwave after it has traveled a distance such that it 
would, on average, hit 1 (solid curve) and 8 (dashed curve) other minihalos, as a function of 
the ambient gas density (depending upon the radii). The two vertical dotted lines enclose 
50% of the minihalos in the large host galaxy. For later discussion (see below), three solid 
dots indicate the velocity of the blastwave from a minihalo with an explosion energy of 
2.6 x 10 54 erg (for e = 0.1) after it has traveled a distance in the mean intergalactic medium 
such that it would have, on average, hit 0.016, 0.13 and 1 (from top to bottom) other 
minihalos, respectively. 
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A higher Mach shock compresses the minihalo gas to a higher density, resulting in 
a lower Jeans mass (roughly inversely proportional to the Mach number). Thus, earlier 
collisions between the faster blastwave and minihalos could lead to formation of lower mass 
stars, even in metal-free gas. We suggest that the observed, extremely low metallicity, low 
mass star HE0107-5240 (Christlieb et al. 2002) may be formed during some early time, 
high Mach shocked minihalos, which are slightly contaminated with metals (see Hideyuki 
& Nomoto 2003). 

One interesting consequence of sweeping blastwave is that a large bulk velocity will 
be imparted on the compressed minihalo clouds. The expected bulk velocity in this case 
is about 3v s /4\^L3 = 0.25t> s , where v s is the blastwave velocity in the ambient gas and 
the factor is due to the density ratio of minihalo gas to the ambient gas. Thus, 

these minihalos may subsequently move at a velocity larger than the escape velocity of the 
large host halo (v esc = 35km s _1 ), if v s > 140km s _1 . Consequently, star formation may be 
widespread spatially and it is even possible that some stars may form after the minihalos 
have left the original host galaxy. One immediate implication is that ionizing photons 
from these stars will be much easier to escape, giving rise to escape fraction possibly well 
in excess of 10%, whereas in the conventional picture of central star formation where the 
escape fraction is estimated to be below 1% (Wood & Loeb 2000; Ricotti & Shull 2001). 

Another implication is that, if some stars have moved outside their original large host 
galaxy when they explode, the resulting blastwaves will travel much further thanks to the 
lower density in the intergalactic medium, which, in turn, could shock heat, accelerate and 
compress minihalos in the intergalactic medium to form stars. If this chain reaction can 
be triggered, star formation in minihalos, which are much more abundant (in total mass) 
than large halos and are otherwise unable to form stars, could produce many more ionizing 
photons than what can be produced by large galaxies alone. In Figure 3 the three solid dots 
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Fig. 4. — shows the evolution of the postshock gas with three different shock velocities, 
v s = (100, 300, 600)km/s, respectively. The panels from top to bottom show the evolution 
of has temperature, density, molecular hydrogen fraction and Jeans mass, respectively. The 
shocks are assumed to propagate through a medium with density equal to the mean gas 
density of the universe at z = 20. 
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show the velocity of the blastwave from a minihalo with an explosion energy of 2.6 x 10 54 erg 
(for e = 0.1) after it has traveled a distance such that it would have, on average, hit 0.016, 
0.13 and 1 other minihalos, respectively (the fraction of mass in minihalos of mass ~ 10 6 M 
fu = 10~ 3 is assumed). If minihalos are strongly clustered with a correlation length 
smaller than 0.2d (where d is the mean separation equal to 14(//10 _3 ) _1 / 3 kpc proper for 
Mm = 10 6 M Q at z — 20), a chain reaction may be sustainable. Minihalos may very well be 
strongly clustered at high redshift (Barkana & Loeb 2003). We will study this separately. 

In any event, blastwaves from stray stars propagating through the intergalactic medium 
would create dense shells. Figure 4 illustrate the temporal evolution of shocked gas at 
z = 20 (assuming isobaric evolution); preshock gas is at the mean gas density at z — 20. We 
see that cooling time is longer in this case than shown in Figure 2; here Compton cooling 
becomes important and its cooling time scale (i.e, e-folding time) is 6 x 10 6 yr. On the longer 
time scales shown in Figure 4 the detailed evolution would be affected by concomitant 
dynamic evolution driven by gravitational instability of large-scale density perturbations. 
Further investigation is needed to see if star formation in shells occurs. It is conceivable 
that a chain reaction may occur, which would result in a picture analogous to the explosion 
scenario proposed for larger galaxies at lower redshift (Ostriker & Cowie 1981). 



4. Conclusions 

In the cold dark matter model large, high-redshift (z > 10) galaxies with virial 
temperature T v > 10 4 K are shown to be mostly comprised of cold atomic clouds, which were 
originally minihalos and accreted into larger halos. It is shown that supersonic collisions of 
cold atomic gas clouds on time scale of 10 7 yr will shock heat and compress the gas, which 
then cools and form H 2 molecules to further cool to ~ 100K. Massive star formation ensues. 
The resulting star formation is spatially distributed. 
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It is then shown that most of the star formation in minihalos may be triggered by 
Shockwaves produced by the first stars formed in collisions. Those shockwave compressed 
minihalos will be more widespread spatially due to large imparted velocities and some can 
possibly escape into the intergalactic medium. If minihalos in the intergalactic space are 
strongly clustered, a chain reaction of star formation in minihalos may be triggered initiated 
by explosions of these stray stars. In this case, far more ionizing photons can be produced 
than large halos alone can. Separately, whether cooling shells of swept-up mean density gas 
in the intergalactic medium can form stars and possibly trigger an explosion-induced chain 
reaction even in the absence of strong minihalo clustering is unclear. 

More importantly with respect to cosmological reionization, widespread star formation 
would allow an ionizing photon escape fraction much higher than in the conventional picture 
where star formation occurs in the central regions of galaxies and whose the escape fraction 
is estimated to be below 1% for galaxies more massive than 10 7 M Q . 

Such star formation may be the norm at z > 10 and may create the physical basis 
to reconcile the standard cosmological model with the WMAP polarization observations. 
An optical depth r e = 0.10 — 0.14 can be achieved realistically even in the absence of 
chain reactions forementioned. A still higher optical depth may indicate that some further 
enhancement in star formation has occurred, such as through possible chain reactions. 
Detailed simulations of star formation triggered by collisions and blastwaves at high redshift 
galaxies would be required to firm up this picture and quantify key parameters, including 
star formation efficiency and ionizing photon escape fraction. 

Finally, we note that these widespread star formation outlined above may give rise to 
metal enrichment of the intergalactic medium far more uniform than any other processes 
that have been proposed (e.g., Gnedin 1998; Aguirre et al. 2001; Cen et al. 2003). 

This research is supported in part by grants AST-0206299 and NAG5-13381. 



-18- 

REFERENCES 

Abel, T., Bryan, G.L., & Norman, M.L., 2002, Science, 295, 93 

Abel, T., Anninos, P., Zhang, Y., & Norman, M.L. 1997, NewA, 2, 181 

Aguirre, A., Hernquist, L., Schaye, J., Weinberg, D.H., Katz, N., & Gardner, J. 2001, ApJ, 
560, 599 

Barkana, R., & Loeb, A. 1999, ApJ, 523, 54 

Barkana, R., & Loeb, A. 2003, astro-ph/0310338 

Bond, J.R., Cole, S., Efstathiou, G., & Kaiser, N. 1991, ApJ, 379, 440 

Bond, J.R., Arnett, W.D., & Carr, B.J. 1984, ApJ, 280, 825 

Bower, R.J. 1991, MNRAS, 248, 332 

Bromm, V., Kudritzki, R.P., & Loeb, A. 2001, ApJ, 552, 464 

Cen, R. 2003a, ApJ, 591, 12 

Cen, R. 2003b, ApJ, 591, L5 

Christlieb, N. et al. 2002, Nature, 419, 904 

Ciardi, B., Ferrara, A., & White, S.D.M. 2003, MNRAS, 344, L7 

Glatzel, W., Fricke, K.J., & El Eid, M.F. 1985, A& A, 149, 413 

Gnedin, N.Y. 1998, MNRAS, 294, 407 

Gunn, J.E., & Gott, J.R., III 1972, ApJ, 176, 1 

Haiman, Z., & Holder, G.P. 2003, ApJ, 595, 1 

Hideyuki, U., & Nomoto, K. 2003, Nature, 422, 871 

Kang, H., Shapiro, PR., Fall, M., k Rees, M.J. 1990, ApJ, 363, 488 

Klein, R.I., McKee, C.F., & Colella, P. 1994, ApJ, 420, 213 



-19- 

Kogut, A. et al. 2003, ApJS, 148, 161 

Lacey, C.G., & Cole, S. 1993, MNRAS, 262, 627 (LC) 

Rakavy, G., Shaviv, G., & Zinamon, Z. 1967, ApJ, 150, 151 

Woosley, S.E. 1986, in Nucleosynthesis and Chemical Evolution, ed. B. Hauck, A. Maeder 
and G. Meynet (Switzerland: Geneva Obs), pi 

Mac Low, M.-M., & Shull, J.M., 1986, ApJ, 302, 585 

Ostriker, J.R, & Cowie, L.L. 1981, ApJ, 243, L127 

Ostriker, J.R, & McKee, C.F. 1988, Rev. Mod. Phys. 60, 1 

Ricotti, M., k Ostriker, J.R 2003, astro-ph/0310331 

Ricotti, M., & Shull, J.M. 2000, 542, 548 

Shapiro, P.R., Iliev, I.T., Raga, A.C., & Martel, H. 2003, in "The Emergence of Cosmic 
Structure", AIP Conference Proc. vol. 666, p. 89 

Smith, J. 1980, ApJ, 238, 842 

Spergel, D.N. 2003, ApJS, 148, 175 

Struck-Marcell, C. 1982a, ApJ, 259, 116 

Struck-Marcell, C. 1982b, ApJ, 259, 127 

Sutherland, R.S., & Dopita, M.A. 1993, ApJS, 88, 253 

Tegmark, M., et al. 2003, astro-ph/0310723 

Wood, K. & Loeb, A. 2000, ApJ, 545, 86 

Wyithe, J.S.B., & Loeb, A. 2003, ApJ, 588, L69 



This manuscript was prepared with the AAS IATgX macros v4.0. 



